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Abstract. The basic model for incompressible two-phase flows with phase 
transitions is derived from basic principles and shown to be thermodynami- 
cally consistent in the sense that the total energy is conserved and the total 
entropy is nondecreasing. The local well-posedness of such problems is proved 
by means of the technique of maximal L p -rcgularity in the case of equal den- 
sities. This way we obtain a local semiflow on a well-defined nonlinear state 
manifold. The equilibria of the system in absence of external forces are iden- 
tified and it is shown that the negative total entropy is a strict Ljapunov 
functional for the system. If a solution does not develop singularities, it is 
proved that it exists globally in time, its orbit is relatively compact, and its 
limit set is nonempty and contained in the set of equilibria. 
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1. Introduction 

Let fl C K™ be a bounded domain of class C 3 ~ , n > 2. ft contains two phases: 
at time t, phase i occupies the subdomain fli(t) of fi, respectively, with i = 1,2. 
We assume that dtti(t) O dil = 0; this means that no boundary contact can occur. 
The closed compact hypersurface T(t) :— dfli(t) C ft forms the interface between 
the phases. 

Let u denote the velocity field, ir the pressure field, T the stress tensor, 9 > 
the (absolute) temperature field, vy the outer normal of fii(t), Vr the normal 
velocity of T(t), Hy = H(T(t)) = — divp^r the sum of the principal curvatures of 
r(i), and [w] = u 2 — v\ the jump of a quantity v across T(t). 
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By the Incompressible Two-Phase Flow with Phase Transition we mean the fol- 
lowing problem: find a family of closed compact hypersurfaces {T(t)}t>o contained 
in SI and appropriately smooth functions u : R+ x fj — > W 1 , ir,9 : M.+ x f2 — > R 
such that 





dtu + u ■ Vu — div T 


= 


in 


n\r(t), 


T-- 


= fi(6) (Vlt + [Vw] T ) - nl, div u 


= 


in 


n\r(t), 


(6){d t e- 


f u ■ V0) - div(<f(0)V0) -T :Vu 
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in 


n\r(t), 




u = d v 9 
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l(0)j + id{6)d v 6] 
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r(t), 




Vt + j -u-vv 
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on 


r(*), 


r(o) 


= r , u(0,x) = u (a), 0(0, a;) 


= (x) 


in 





(1.1) 



The variable j, called phase flux, can in fact be eliminated from the system: solving 
for j yields 

j = -id{p)d v oyi{o), 

provided 1(9) ^ 0, a property which is assumed later on anyway. In fact, if 1(0) = 
then the interfacial velocity V is not uniquely determined which leads to non-well- 
posedness of the problem. The equation for Vr then becomes 

Vr-u -vr- \d(e)d v eyi(e) = o. 

Note that the sign of the curvature Hr is negative at a point x G T if Oi n B r (x) 
is convex, for some sufficiently small r > 0. Thus if Cl% is a ball, i.e. T = Sr(xq), 
then H r = -(n - 1)/R. 

Several quantities are derived from the specific Helmholtz free energies ipi (&) as 
follows: 

• Ci(6) = i/ji(6) + Grji(6) means the specific internal energy in phase i, 

• Vi(0) = — ^(0) the specific entropy, 

• Ki{6) = e-(6*) > the specific heat capacity, and 

• 1(9) = e\4,'(6)] = -0[r?(0)] the latent heat. 

Further di(9) > denotes the coefficient of heat conduction in Fourier's law, 
fJ,i(9) > the viscosity in Newton's law, p :— p\ = p2 = 1 the constant density, 
and a > the constant coefficient of surface tension. In the sequel we drop the 
index i, as there is no danger of confusion; we just keep in mind that the coefficients 
depend on the phases. 

This model is explained in more detail in the next section. It has been recently 
proposed by Anderson et al. [2], see also the monographs by Ishii [5] and Ishii 
and Takashi [10]. We will see below that it is thermodynamically consistent in 
the sense that in the absence of exterior forces and of fluxes through the outer 
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boundary, the total energy is preserved and the total entropy is nondecreasing. It 
is in some sense the simplest sharp interface model for incompressible Newtonian 
two-phase flows taking into account phase transitions driven by temperature. 

There is an extensive literature on isothermal incompressible Newtonian two- 
phase flows without phase transitions, and also on the two-phase Stefan problem 
with surface tension modeling temperature driven phase transitions. On the other 
hand, mathematical work on two-phase flow problems including phase transitions 
is rare. In this direction, we only know of the papers by Hoffmann and Starovoitov 
IB], dealing with a simplified two-phase flow model, and Kusaka and Tani [T31 
If 4] which are two-phase for the temperature but with only one phase moving. 
The papers of DiBenedetto and Friedman [3] and DiBenedetto and 0'Leary[4] 
deal with weak solutions of conduction-convection problems with phase change. 
However, the models considered in these papers do not seem to be consistent with 
thermodynamics. 

It is the purpose of this paper to present a rigorous analysis of problem (jl.ip in 
the framework of L p -theory. We prove local well-posedness, and we identify the 
equilibria of the problem. In turns out that the equilibria are the same as those for 
the thermodynamically consistent two-phase Stefan problem with surface tension. 
This heavily depends on the fact that the densities of the two phases are assumed 
to be equal; in this case the problem is temperature dominated. 

In a forthcoming paper we will consider the case where the densities are not 
equal; then the solution behavior is different, as the interfacial mass flux has a 
direct impact on the velocity field of the fluid. Physically this relates to the so- 
called Stefan currents which are induced by phase transitions. The velocity field is 
then no longer continuous across the interface and this leads to different analytic 
properties of the model. We call this case velocity dominated. 

The plan for this paper is as follows. In the next section we introduce and dis- 
cuss the model in some detail and show that it is thermodynamically consistent. In 
Section 3 we prove that the negative total entropy is a strict Ljapunov functional 
for the problem, i.e. is strictly decreasing along nonconstant solutions. We iden- 
tify the equilibria which consist of velocity zero, constant temperature, constant 
pressures in the phases, and the disperse phase Oi consists of finitely many nonin- 
tersecting balls of the same radius. This is in contrast to the isothermal two-phase 
Navier-Stokes equation with surface tension where phase transitions are neglected: 
in that case the balls may have different sizes! 

In our approach we apply the well-known direct mapping method where the 
problem with moving interface is transformed to a problem on a fixed domain. This 
transformation will be introduced in Section 4; as a result we obtain a quasilinear 
parabolic evolution problem with dynamic boundary condition. In Section 5 we 
study maximal L p -regularity of the underlying linearized problem. It turns out 
that the linear problem is equivalent to the linearized Stefan problem with surface 
tension, followed by a Stokes problem where the curvature serves as an input 
variable on the interface. Here we rely on previous work Priiss, Simonett and 
Zacher [21]. Local well-posedness is discussed in Section 6, while the proof of the 
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main result is given in Section 7. It is based on maximal L p -regularity of the linear 
problem and the contraction mapping principle. Finally, in Section 8 we consider 
global existence of solutions. It is shown that solutions exist globally in time, 
provided the interface satisfies a uniform ball condition, the temperature and the 
modulus of the latent heat are bounded from below, and certain a priori estimates 
on the solution hold. In a forthcoming paper we will show that in this situation 
the solution converges to an equilibrium in the topology of the state manifold. 

2. The Model 
In this section we explain the model. We begin with 
Balance of Mass: 

d t p + div (pit) = in 0\T(i), 

[p(u - u r )] • W = on T(t), 

where ur is the interfacial velocity. Hence the normal velocity Vr can be expressed 
by Vr = Ur ■ vv- We define the phase flux, more precisely, the interfacial mass 
flux, by means of 

j := p(u - u r ) ■ vr, i.e. = [u • vrj, (2.1) 

and we note that j is independent of the phases, hence well-defined. Therefore, a 
phase transition takes place if j ^ 0. If j = then u ■ Vr = Ur • vr = Vr and in 
this case the interface is advected with the velocity field u. 

In this paper we consider the completely incompressible case, i.e. we assume that 
the densities are constant in the phases fij. Then conservation of mass reduces to 

divu^O in Q\T(t). (2.2) 

If only the latter property holds, we say that the material is incompressible. 

Integrating over the finite domain Sli(t) we obtain by the Reynolds transport 
theorem 

^- ( pdx = [ pV~rds+ ( d t pdx 
dt Jnut) Jrtt) Jf2i(t) 



/ pVr ds — div(pu) dx 

Jr(t) Jntit) 

/ (pur ■ vr — pu ■ vr) ds = — I j ds. 
Jrtt) Jr(t) 



T(t) Jr{t) 
In case |^ 2 (£)l i s finite as well, we obtain in the same way 

d_ 
~dt 



pdx = j j ds, 
n 2 (t) Jr(t) 



proving conservation of total mass. 

Therefore, if both phases are completely incompressible we have 

Pi|fti(t)| + P2|n 2 (t)| = Pi|«i(0)| + p 2 |O 2 (0)| =: c . 
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This implies 

[p]|fii(t)| = P2|n|-co, 

hence |f2i(t)| = constant in the case of nonequal densities, i.e. the phase volumes 
are preserved. If the densities are equal, there is no restriction on the phase 
volumes due to conservation of mass. This constitutes a big difference bctwecen 
the two cases. 

In later sections of this paper we are concerned with the completely incom- 
pressible case where, in addition, the densities are equal, and hence conservation 
of mass does then not imply conservation of phase volumes. 

Next we consider 

Balance of Momentum: 

d t (pu) + div (pu <g> u) -divT = pf in 

\pu ® (u - u r ) - T\v T = div r T r on T(t). 

Using balance of mass and the definition of the phase flux j we may rewrite this 
conservation law as follows: 

p(d t u + u- Vit) -divT = pf in Q\T(t), 

Mi - [Tut] = divrTr on T(t). [2 " i} 

Here Tr denotes surface stress, which in this paper will be assumed to be of the 
form 

T r = erP r , Pt = I-v t ®v t , 
where a > is a constant. So in particular we neglect here surface viscosities as 
well as Marangoni forces. Then we have 

div r T r = oH v v v (2.4) 
which is the usual assumption for surface stress. 
We now turn our attention to 

Balance of Energy: 

Let q denote the heat flux and r an external heat source. Then balance of energy 
becomes 

d t (^\u\ 2 2 +pe) + div{(^\u\ 2 2 + pe)u}-div(Tu-q)=pf-u + pr in 

{{^\u\l+pe){u-u T )-{Tu-q)\-v T = div r T r ■ u T on T(t). 

Note that the second line implies that the only surface energy taken into account 
here is that induced by surface tension. 

Using balance of mass, balance of momentum, and the definition of the phase 
flux j we may rewrite this conservation law as follows: 

p(dte + u ■ Ve) + div q — T : Vu = pr in 
(lej + ll\u-u r \l])j-lTv r .(u-u r )] + lq-vr]=0 on T(t). 
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The total energy is given by 

E:=E(u,9,T) :=\ I p\u\\dx + [ pe{9) dx + a\T\, 
2 Jn Jn 

where |T| denotes the surface area of V. Note that here the assumption of a being 
constant is important! Otherwise, one has to take into account surface energy, 
locally, which means that we would have to include a balance of surface energy. 

In the absence of external forces / and heat sources r, for the time derivative 
of E we obtain by the transport theorem 

j t E = j {u ■ pd t u + pd t e(6)} dx - J{[^\u\\ + pe{6)\ + aH r }V r ds 

{p(u ■ V)m • u — div T ■ u + p(u ■ V)e(0) + div q — T: Vw} dx 

~ J {[f Hi + pe{0)] + aH T }u T ■ v v ds 

= J {[(f M2 + P<0))(u - u T ) ■ M - [Tv T ■ uj + [q ■ ur] - *H r V r } ds 

= f r {l\\ u \ 2 2 + <°)y - l Tu ■ M + ll ■ M - CTff rVr} ds, 
provided that energy transport through the outer boundary is zero, which means 
old 2 

{i-^-+pe{9)}u-u-Tvu + q-v = on dVl. 

Hence 

!E(M,r) = o, 

which implies that the total energy is preserved. 
Entropy Production: 

We introduce now the fundamental thermodynamic relations which read 

e(0) = v(0) + <W), v (d) = -^(0), 

k(6) = e'{6), 1(6) - 9[1/(0)] = -0[v(6)], 

where ip(9) means the Helmholtz free energy which should be considered as given, 
but depends on the phases. The quantities r\ and n, I are called entropy, heat 
capacity, and latent heat, respectively. 
The total entropy is defined by 

$(0,T) = / pn(0)dx. 

The main idea of the modeling approach presented here is that there should be no 
entropy production on the interface, as it is considered to be ideal; in particular 
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it is assumed to carry no mass and no energy except for surface tension. We have 
with (|23 ]) - (|2l3]l for / = r = 0, and with e'{9) = <V(6>), 



*(u,0,r)= f pd t r,{9)dx- [ Ip^jVrds 
Jn Jr 

pv'(e)d t e dx- j \pn(0)]u T ■ w ds 



d_ 

Tit 



= L { wi {T : D ~ div q} ~ pu ' Vv{e) } dx ~ I Ip7?(0)1 ur ' vt ds 

= f u {^r- ^} dx + J r {lPV(°)(u u r ) ■ u r ] + l[q ■ ^r]} da 

= /"{^-i^}^ + /J{[- W + b 

provided there is no entropy flux through the outer boundary dQ, which means 

q ■ v + p6r)(8)u ■ v = on 9f2. 

To ensure conservation of energy as well as entropy through the outer boundary 
we impose in this paper for simplicity the following 

Constitutive Laws on the Outer Boundary dfl: 

q-v = 0, u = 0. (2.6) 

Actually, we may consider more general conditions at the outer boundary, like 
q ■ v = and the partial slip condition u ■ v = 0, PqqTv = 0, but we refrain from 
doing this here. 

As constitutive laws in the phases we employ Newton's law for the stress tensor 
and Fourier's law for the heat flux; these ensure nonnegative entropy production 
in the bulk. Recall that we consider the completely incompressible case. 

Constitutive Laws in the Phases: 

T = S-irI, S = 2fi(9)D, D = i(Vu + [V U ] T ), 
q = -d{6)V9. 

We assume k(8) = e'(9) = —9ip"(9) > 0, as well p(9) > for the shear viscosity, 
and d{9) > for the heat diffusivity; note that each of these quantities may jump 
across the interface. 

Finally, on the interface we assume 

Constitutive Laws on the Interface T(t): 

[0] = [P r u]=O, -l(9)j + \q ■ v T \ = 0. (2.8) 

Hence in this model the temperature and the tangential part of the velocity are 
continuous across the interface; the latter means that there is no tangential slip 
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at the interface, as it is considered ideal. The third equation in (|2.8|) means that 
no entropy is generated on the interface. This implies the relation 



dt 



mr)= [ {^\D\l + ^e\l} dx >Q (2.9) 

for the entropy production. 

Employing the constitutive assumption \Pvu\ = with ftpVr] = yields 
P 

{Tit ■ (u - ut)] = [Tur ■ Pr(u - u r )j + {-Tit ■ vr\j = \-Tv T ■ v T \j. 

Combining these conditions with the energy balance across the interface and the 
constitutive law for q results into the following jump conditions, assuming that j 
may take arbitrary values: 

l(6)j + {d{6)d v 6\ = on T{t), 



l,o,l T1 , , (2.10) 



{W)] + Uj 2 -Kfr-H}j=0 on r(t). 



The first of these two equations is the Stefan law and the second the (generalized) 
Gibbs- Thomson law. 

Note that if j = 0, i.e. in the absence of a phase transition, the Stefan law 
becomes [d(0)d„0] = 0, and the Gibbs- Thomson relation trivializes. On the other 
hand, if phase transitions occur then j must be considered as arbitrary, and so the 
Gibbs- Thomson law becomes 

W)\ + 1 Ab' 2 - i-Tvv ■ M = o- (2-11) 

2p z p 

We summarize to obtain the following model for incompressible two-phase flows 
with phase transitions. 

The Complete Model 

In the bulkQ\T(t): 

p(d t u + u-Vu)-divT=0, 
T = /i(0)(Vu + [Vu] T ) — ttI, divu = 0, (2.12) 
pn(9)(d t 6 + u ■ V0) - div(d(0)V0) — T : Vu = 0. 



On the interface T(t): 

,1, 



Itt^Jj 2 - l-Tir -»rl = = "• » i : " ' "j J- 

2p P P 



0. 






M 


= \-\jvr 






P 


0, 


m 


= 0, 


0, 




= u ■ IT - 



(2.13) 
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On the outer boundary dQ: 



u = 0, d v 9 = 0. 



(2.14) 



Initial conditions: 



r(o) = r , «(o) = mo, 0(0) = e Q . 



(2.15) 



This results in the model (jl.ip if we assume in addition that the densities are 
equal, i.e. pi = pi = 1. 



As we have seen in Section 2, the negative total entropy is a Ljapunov functional 
for the problem, and it is even a strict one. To see this, assume that $ is constant 
on some interval (h,t 2 ). Then j- t <&(9, V) = in (ti,t 2 ), hence D = and V0 = 
in (ti,ti) x fi- Therefore, 9 is constant which implies \d(9)d v 9\ = 0, and then 
from the interfacial boundary condition we obtain j = 0, provided 1(9) ^ 0. Then 
\u\ = 0, hence by Korn's inequality we have Vu = and then u = by the no-slip 
condition on dfl. This implies further dtO = dtu — and Vr = 0, i.e. we are 
at equilibrium. Further, V7r = 0, i.e. the pressure is constant in the components 
of the phases, and aHr — fnj, {ip(6)j = — W are constant as well. Since 9 is 
continuous across the interface the last relation shows that it is constant in all of 
fil, even if it is not connected. From this we finally deduce that fii is a ball if 
it is connected, or a finite union of balls of equal radii, as fi is bounded. Let us 
summarize. 

Theorem 3.1. Let a > 0, £ C 3 (0, oo), pi, di € C 2 (0, oo), 

-^(s),p l ( S ),d l ( S )>0, [V'(s)]^0, 

for s > 0, i = 1, 2. TTien i/ie following assertions hold. 

(i) TTie <o<aZ energy E is conserved along smooth solutions. 

(ii) TTie negative total entropy — $ is a strict Ljapunov functional, which means 
that — $ is strictly decreasing along nonconstant smooth solutions. 

(iii) TTie non- degenerate equilibria (i.e. no boundary contact) are zero velocity, 
constant temperature, constant pressure in each phase, and f2i consists of 
a finite number of nonintersecting balls of equal size. 

Therefore the equilibria are the same as those for the thermodynamically con- 
sistent Stefan problem with surface tension which has been discussed recently in 
Priiss, Simonett and Zacher |21) . We now discuss some basic considerations from 
that paper. Suppose (u, 7r, 9, V) is an equilibrium with 



3. Equilibria 
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where B(x, R) denotes the ball with center x and radius R, and S(x, R) its bound- 
ary. To determine ([7r], 0, R) we have to solve the system 

<p(0,R) := \n 1 \e 1 (9) + \n 2 \e 2 (9)+a\T\=E , 

W = aH r , (3.1) 

[#?)] = -M, 

where Eo = E(uo,6o,Ff)) is the prescribed initial total energy, and <p(9,R) means 
the total energy at a given equilibrium to be studied. We have 

[tt] = aH r = -(n - l)a/R, R = R(9) = (n - l)<r/ty(0)]. 

Hence there remains a single equation for the equilibrium temperature 9, namely 

E = <p{6) := Me 2 (9) - m(u Jn /n)R n (9){e{9)\ + amtu n R n - 1 (9), 

where u) n denotes the area of the unit sphere in R™. Note that only the temperature 
range [?/>(#)] > is relevant due to the requirement R > 0, and with 

R* m = sup{i? > : VL contains m disjoint balls of radius R} 

we must also have R < R^, i-e. with h{9) = \ip{9)\ 

m > 

With e(9) = ip(0) — 9ip f (Q) we may rewrite (p(9) as 

^)H% 2 W + c.(^ + (n-l)^), 

where we have set c„ = m— — - — -((n — l)a) n . 

n(n — 1) 

Next, with 

a{n-l)h'{6) h'(9)R 2 {9) 



R'{9) 



h 2 {9) cr(n-l) 



we obtain 

<p'{9) = \Q\e' 2 (9) [e'(0)]|fii| + mw „(^--il - [e^)])^- 1 ^)^) 



'v{n- 1) 

\n\n 2 (9) - [K(0)]|fti| +muj n 9h'{9)R n - 1 {9)R'{9) 
.h'{9)R 2 {9) 



= {K\\) a -9h'{9)\T\- 



a(n - 1) 

(k|1)oJ? 2 (0) f a(n-l) Z(6>) 2 |r| ^ 
a(n-l) I R 2 {0) 0(«|l) n J' 



with = 9h'(9). It will turn out that the term in the parentheses determines 
whether an equilibrium is stable: it is stable if and only if m = 1 and <J>'(9) < 0. In 
particular, if fli is not connected this equilibrium is unstable; this fact is related 
to the physical phenomenon called Ostwald ripening. 
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In general it is not a simple task to analyze the equation for the temperature 

unless more properties of the functions £2(0) and, in particular, of [V'(^)] are 
known; cf. Priiss, Simonett and Zacher |21j for further discussion and results. 

4. Transformation to a Fixed Domain 

Let f2 C R n be a bounded domain with boundary d£l of class C 2 , and suppose 
r C £1 is a closed hypersurface of class C 2 , i.e. a C 2 -manifold which is the boundary 
of a bounded domain fix C f2; we then set = Note that H2 is connected, 

but Hi maybe disconnected, however, it consists of finitely many components only, 
since dili = T by assumption is a manifold, at least of class C 2 . Recall that the 
second order bundle of T is given by 

A/" 2 T := {(p,ur(p)>^rMp)) '■ P e F}. 

Here vr{p) denotes the outer normal of Vt\ at p £ T and Vr the surface gradient 
on T. The Wcingarten tensor Ly on T is defined by 

Lr(p) = -Vr^r(p), p £ T. 

The eigenvalues Kj(p) of Lr(p) are the principal curvatures of T at p £ T, and we 
have |ir(p)| = max.,- The curvature Hr(p) (more precisely (n — 1) times 

mean curvature) is dehned as the trace of Lr(p), i.e. 

n-i 

Hrip) = ^2 K i( p ^ = - div r^r(p), 
j'=i 

where divr means surface divergence. Recall also the Hausdorff distance dn be- 
tween the two closed subsets A, B c R m , defined by 

dif(A,B) := maxjsup dist(a, -B), sup dist(6, A)}. 

a£A beB 

Then we may approximate T by a real analytic hypersurface E (or merely E G C 3 ), 
in the sense that the Hausdorff distance of the second order bundles of T and E is 
as small as we want. More precisely, for each r\ > there is a real analytic closed 
hypersurface E such that d//(A/" 2 E,7V 2 r) < n. If i] > is small enough, then E 
bounds a domain fif with fif C fi, and we set J7f = . 

It is well known that a hypersurface E of class C 2 admits a tubular neighbor- 
hood, which means that there is a > such that the map 

A : E x (-a, a) -> M" 
A(p,r) := p + ruY,{p) 
is a diffcomorphism from E x (—a, a) onto 7?.(A). The inverse 

A -1 : 11(A) h4 E x (-a, a) 
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of this map is conveniently decomposed as 

A- 1 (.x) = (n s (x),d s (x)), xeK(A). 

Here lis (a;) means the orthogonal projection of x to £ and <fe(x) the signed 
distance from x to E; so c?s(x)| = dist(x, E) and ds(x) < if and only if a; € Clf . 
In particular we have 11(A) = {x G W l : dist(x, E) < a}. 

Note that on the one hand, a is determined by the curvatures of E, i.e. we must 
have 

< a < min{ \/\nj(p)\ : j = 1 , . . . , n — 1, p G £}, 

where «j(p) mean the principal curvatures of £ at p G £. But on the other hand, 
a is also connected to the topology of £, which can be expressed as follows. Since 
£ is a compact C 2 manifold of dimension n — 1 it satisfies the ball condition, which 
means that there is a radius r E > such that for each point p G £ there are points 
Xj G Of, i = 1,2, such that B rs (xj) C fif , and B rs (xi) fl £ = {p}. Choosing r E 
maximal, we then must also have a < r E . In the sequel we fix 

a = imin{r s , 1 : j = 1, . . . , n - l,p G £}. 
2 |Kj(p)| 

For later use we note that the derivatives of lis (x) and c?s (x) are given by 

W s (x) = v£(n E (aO), Dn E (s) = M (cfe(z))(n(x))ft(n E (x)) 

for |d E (x)| < a, where P E (p) = I — vy,(p)®vi:(p) denotes the orthogonal projection 
onto the tangent space T p Y> of £ at p G £, and 

M (r)(p) = (J-rLstp))- 1 , (r,p) G (-a, a) x £. 

Note that 

|M (r)(p)| < 1/(1 - r|L E (p)|) < 2, for all (r,p) G (-a, a) x £. 

Setting r = T(t), we may use the map A to parameterize the unknown free bound- 
ary T(t) over E by means of a height function h(t,p) via 

r(t) = {p + ft(t,pK(p): P es, t>o}, 

at least for small \h\oo. Extend this diffeomorphism to all of f2 by means of 

E h (t,x) = x + x(ds(x)/a)h(t,n s (x))u T ,(Il T ,(x)) =: x + £ h (t,x). 

Here \ denotes a suitable cut-off function. More precisely, \ G V(M.), < \ < 1, 
X(r) = 1 for |r| < 1/3, and xW = for \r\ > 2/3. Note that E h (t,x) = x for 
\d(x)\ > 2a/3, and 

x) = x — h(t, x)vy,(x), x G E, 
for \h\oo sufficiently small. Now we define the transformed quantities 

u(t, x) = u(t, £/j(t, a;)), n(t,x) — Tr(t,Eh(t,x)), t > 0, x G fi\E, 
0(i, x) = 0(i, S h (i, x)), f > 0, x G 0\E, 
j(t,x) =j(t,E h (t,x)), t>0, xgE, 
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TZ U (u,6,h) 


in 







in 


n\s, 


Ke(u,9, h) 


in 


n\s, 





on 


asi, 


aH v {h)v T (h) 


on 


s, 





on 


s, 





on 


s, 





on 


s, 


/lo- 







the ptdZ feacfc of (u,tt,9) and j. This way we have transformed the time varying 
regions \ r(i) to the fixed domain \ S. 

This transformation gives the following problem for (u, n, 8, h): 

d t u - Q{h) ■ »(6)(g(h)u + [g(h)u] T ) + g(h)Tt 

G{h) ■ u 

K(0)dt6 - G(h) ■ d(6)G(h)6 
u = d v Q 

-lv{0){G{h)u+[G(h)u] T )Mh)-nMh)l = <jH T {h)v T {h) on S, (4.1) 

\u\ = n 

[iP(e)}+aH r (h) 

/3(h)d t h -U-W- {d{e)G(h)e ■ M/m 

u(0) = u , 0(0) = S , h(0) 

Here G(h) and Hr(h) denote the transformed gradient and curvature, respectively. 
More precisely we have 

D~ h = I + D£ h , [D E^- 1 =/-[/+ D^D^ =:I- M^hf , 

with 

D£ h (t,x) = ^x'(d^(x)/a)h(t,U^(x))^(U^(x)) «>^ E (n E (x)) 
+ X(d s (x)/a)^ s (n s (x)) ®M (d E (a;))V E ^(t,n s (x)) 
- X (d E {x)/a)h(t, IT E (x))L E (n(x))M (d E (x))P E (n E (x)) 
-/ l (t,n(x)) J L s (n(x))M (d s (x))P s (n s (x)) for |d E (aO| < 2a/3, 

D£ h (t,x) = for |d s (a) | > 2a/3. 
In particular, 

D£ h {t,x) = ^s(n(x)) <8»Mo(d E (a;))V E /i(t,n(a;)) for |d E (x)| < a/3. 
Thus, [/ + D£h] is boundcdly invertible if h and V E ft. arc sufficiently small, e.g. if 

Noo< ^min{a/| X '|oo,l/|£s|oo} and (Vs/i^ < 1. (4.2) 
With these properties we derive 

VioH^ 0(fc)7f = [VS^ 1 o E h ] T V7r = [VS^-^Vtt = (/ - M 1 (/i))Vtt 
divu o S/j = C/(fc) • u = (I - Mi(/i))V • tt. 
Next we note that 

d t u oE h = d t u - [Du o E h ]d t E h = d t u - Du [D S^]" 1 d t £h 

= d t u - {[I + Dth}- 1 dt£ h • V)u =: d t u + (R(h) • V)«. 
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with R(h) = -[1+ D^h}- 1 dtth- Hence 

TZ u (u, 8, h) = -u- 9(h)u - (R(h) ■ V)u. 

Observe that the function R(h) contains a time derivative of h linearly. Similarly 
we get 

V0 o E h = Q(h)6 = (I- Mi(/i))V0 

d t 9 o E h = 9t§ - V0 ■ [VE^dtth =: d t 9 + (R(h) • V)0, 

and so 

TZ e (u, S,h)) = -k{6)u ■ Q(h)9 - k(6) (R(h) ■ V)0 + fi(9) (Q(h)u + [G(h)uf) : Q(h)u. 
With the Weingarten tensor Ls and the surface gradient Vs we further have 
Mh) = PWfa - a(h)), a(h) = M (h)\7 s h, 
M (h) = (7- hL^)-\ (3(h) = (1 + \a(h)\ 2 )- x '\ 

and 

Vr = d t E ■ v v = d t hv r ■ = /3(h)d t h. 

The curvature Hr{h) becomes 

H T {h) = /3(/i){tr[M (/i)(L E + V s a(/i))] - (3 2 (h)M (h)a(h) • [V s a(/i)]a(/i)}, 

a differential expression involving second order derivatives of h only linearly. Its 
linearization is given by 

/ff(0)=tri| + A s . 

Here denotes the Laplace-Beltrami operator on S. 

It is convenient to decompose the stress boundary condition into tangential and 
normal parts. For this purpose let as before Py, = I — v-^^vy, denote the projection 
onto the tangent space of X. Multiplying the stress interface condition with vt,/ ft 
we obtain 

[if] - aH r (h) = M6)(g(h)u + [g(h)u] T )(^ - Mo(h)V^h) ■ i^j, 

for the normal part of the stress interface condition. Substituting this expression 
for [fr] — aHr(h) in (|4.1[) and applying the tangential projection yields 

- PzM8)(g(h)u + [g(h)u] T )(^ - M {h)Vvh)j 

= W){${h)u + [G(h)u] T )(w - M (h)Vvh) ■ usWoWVvh 

for the tangential part. Note that the latter neither contains the pressure jump 
nor the curvature which is the advantage of this decomposition. 

The idea of our approach can be described as follows. We consider the trans- 
formed problem (|4.1|L Based on maximal L p -regularity of the linear problem given 
by the left hand side of (14. ip , we employ the contraction mapping principle to ob- 
tain local well-posedness of the nonlinear problem. This program will be carried 
out in the next sections. 
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5. The Linear Problem 
The principal part of the linearization of (|4.1j) reads as follows. 

dtu — fj,o(x)Au + V7r = f u in fi\E, 

div u = fd in 57 \ E, 

[uj = on E, 

-lpo(x)(Vu + [Vu] T )i/s - 7ri/ s ] = aA^hvy, + g u on E, 

u = on c>f2, 

u(0) = uq in f2, 

K (x)d t 6 - d (x)A6> = /„ in ft \ E, 

[6>] = on E, 

9„6> = on 90, 

6>(0) = 6» in ft, 



(5.1) 



(5.2) 



Here 



h(t, x)6 + (tAj/i = ge on E, 

9th-[do(x)du0\Mx)=9h on E, (5.3) 
/i(0) = h on E. 

// (a;) = fi(6 (x)), k (x) = k[9 [x)), d (x) = d(6 (x)), 

l (x) = l(6 (x)), h(t, x) = [^((e Ast 7^o)(^))], 

where jOq means the restriction of 9q to E. Observe that the term u ■ vy in the 
equation for h is of lower order as it enjoys more regularity than the trace of 6 on 
E. Therefore this system is triangular, since u does neither appear in (|5.2| nor 
in (|5.3p . The latter system comprises the linearized Stefan problem with surface 
tension which has been studied in [3TJ Theorem 3.3]. To state this result we define 
the solution spaces 

E e (J) = {6 e Hp(J: L p (Q)) n L p ( J; H%(Q \ E) n C(ft)) : ^6» = on 90}, 

MJ) = W^'^iJ; L P (E)) n W p 1 ~ 1 l 2p (J\ tf p 2 (E)) n L p (J; % 4 -^(E)), 

E(J) =E 9 (J) xE ft (J) 

and the spaces of data 

F fl (J) = L p (J x 0), 

= W^- 1/2p (J; £„(£)) n i p (J; % 2 - 1/p (E)), 
¥ h (J) = wV2-i/2p ( j. Lp(E) ) n L P (J; Wp 1 " 1 ^ (S)), 
F(J)=F 9 (J)xF H (J)xF ft (J). 

Then we have 
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Theorem 5.1. Let p > n + 2 and a > 0, suppose kq € C(Cli), do £ C 1 ^), 
i = 1,2, Ko,do > on Cl, Iq G C^E), and let l\ G ¥h(J) such that IqI\ > on 
JxS, where J = [0,to] *s a finite time interval. 

Then there is a unique solution (9, h) G E( J) o/ (|5 . 2|) - (|5 .3[) if and only if the data 
(fe,ge,9h,0 Q ,h Q ) satisfy 

(fe,9e,9h)^HJ), (0o, h Q ) G [W p 2 -^(fl\S)nC(ii)] x W^(Z), 
and the compatibility conditions 

d v = on dn, Zi(O)0 o + vAxho = g e (0), g h {0) + [doft/^ol/Jo 6 W^ 6 / p (£). 

T7ie solution map [(fo, go, gh, Za) l ~ ?• (^i^ 1 )] ^ continuous between the corresponding 
spaces. 

Having solved this problem, Ag/i G J), hence we may now solve the remain- 
ing Stokes problem to obtain a unique solution (u, 7r) in the class 

u G flJ(J; L p {tt)) n n L p (J; H 2 p {VL \ s) n C(Q)) n , Vtt G L P (J, L p (n)) n , 

provided the data satisfy 

f u e¥ e (J) n , g u e¥ h (J) n , /^^(J^^nL^^f^s)), 

as well as uq G \Wp~ 2 ^ p (Cl \ S) fl C(f2)] n , and the compatibility conditions 
u = 0ondil, dwu o = f d {0), -Psbo(Viio + [Vii ] T )^] =Ps.9«(0), 

are satisfied. Concerning the two-phase Stokes problem we refer to Shibata and 
Shimizu [22] and to Kohne, Priiss and Wilke [IT], where the case fio = const 
has been treated. This result can be extended to nonconstant fiQ via the method 
of localization. For the one-phase case this has been carried out for much more 
general Stokes problems in Bothe and Priiss [JJ 

Therefore the linearized problem (|5.1I) - (15.3[) has the property of maximal Ir- 
regularity. To state this result we set 

E„(J) :={ue [H^(J;L p (n))r\L p (J;H^(n\T,)r\C(Q))] n : u = on dQ}, 
E n (J) :=L p (J;flJ(n\E)), 

and define the solution space for (|5.1|) - (|5.3j) as 

E(J) = E U (J) x E W (J) x Fft(J) x Ee(J) x E h (J). 

Then the main result on the linearized problem reads 

Theorem 5.2. Let p > n + 2 and a > 0. Suppose ^o, € C(f^)> do G C 1 (^i) 7 
i = 1,2, kq, do > on f2, G C (£<), h G Eft(J) swc/i £/iai Zo^i > on J xY,, 
where J = [0, Iq] *s h finite time interval. 
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Then the linear problem (I5.ip - (I5.3|) admits a unique solution (u,tt, [7r],#, h) G 
E(J) if and only if the data (uo,0o,ho) and (f u , fd, 9u, fe,ge,9h) satisfy the regu- 
larity conditions: 

(u , ho) G % 2 - 2/p (fi \ E)» x w£- 2 /p (fl \ S) x <- 3 /p(E), 

(A,, /«) e F e (j) n+1 , / d e Bjfj; H-\n)) n l p ( j ; ff^n \ E)), 
(9u,g h )e¥ h (J) n+1 , g e e¥ H (j), 

and the compatibility conditions: 

divu = fd(0) in Q\E, 

u = d u 9 = on dfl, 

M = M=0 on E, 

-P s [^o(Vwo + [Viio] T )^] =,Ps.9„(0) on E, 

Zi(O)0 o + <TA s /i O =5e(0) on E 
0,(0) + [doMo]/io G % 2 - 6 / p (E). 

TTie solution map [(u ,8 ,ho, f u , fd, 9u, fe, 99, 9h) >-> (ti,!r,[4M)] continuous 
between the corresponding spaces. 

6. Local Well-Posedness 

The basic result for local well-posedness of problem in an L p -setting is the 
following. 

Theorem 6.1. Let p > n + 2, a > 0, suppose ipi G C 3 (0,oo) 7 (ii,di € C 2 (0, oo) 
such that 

m(s) — -stp"(s) > 0, (ii(s) > 0, dj(s) > 0, s G (0, oo), i = 1, 2. 

Lei fi C K n be a bounded domain with boundary <9fi G C 3 ~ and suppose To C f2 
is a closed hyper surf ace. Assume the regularity conditions 

u g w^ 2 /p(n \ r ) n , e g wl~ 2 'v (fi \ T ), To g w^, 

the compatibility conditions 

divwo = in fi \ To, 

no = d v 9o = on <9fi, 

M = [0 O ]=O on r„, 

Fro[Mo(VMo + [Vuo] T )tTo] =0 on T , 

Meo)]+<rH To =0 on T , 

[doMo] g w^- 6/p (r ), 

and i/ie well-posedness condition 

2(#o) 7^ on Tq and Oo > on fi. 
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Then there exists a unique L p -solution of problem (jl.ljl on some possibly small but 
nontrivial time interval J = [0,a]. 

Here the notation Tq 6 Wp 3 ^ p means that Tq is a C 2 -manifold, such that its 
(outer) normal field vr is of class Wp~ 3 ^ p (r Q ). Therefore the Weingarten tensor 
L Ta = -V ro ^r of T belongs to Wp~ 3/p (T ) which embeds into C 1+a (r ), with 
a = 1 — [n + 2)/p > since p > n + 2 by assumption. For the same reason we also 
have u G C 1+Q (^(0)))", and 9 G C 1+Q (^(0))), i = 1,2, and V G C 2Q (r ). 
The notion L p -solution means that (u, 7T, 0, T) is obtained as the push-forward of 
an .Lp-solution of the transformed problem (14. The proof of Theorem 16.11 is 
given in the next section. 

For later use we discuss an extension of the local existence results to spaces with 
time weights. For this purpose, given a C/MD-Banach space Y and /i G (1/p, 1], 
we define for J = (0, to) 

K S P JJ;Y) := {u G L Ptloc (J;Y) : ^""u G ^(J;^)}, 

where s > and if G {if, 14^}. It has been shown in Priiss and Simonett [17] that 
the operator d/dt in L Ptfl (J;Y) with domain 

D(d/dt) = oHlJJ; Y) = {ue H^J; Y) : u(0) = 0} 

is sectorial and admits an _ff°°-calculus with angle tt/2. This is the main tool to 
extend the results for the linear problem, i.e. Theorem I5.2| to the time weighted 
setting, where the solution space E(J) is replaced by E^(J) and F(J) by F M (J) 
respectively, where 

z G E M (J) O G E(J), / G F M (J) <S> t x -^j G F(J). 

The trace spaces for (u, 6) and h for p > 3 are then given by 

(u ,6» ) G [Wp^ 2/p {Q\T.)r\C(Cl)] n+1 , fcoeff^lE), 

fti := (ftft)(0) G % 4 ^ 2 - 6 /P(S), !(> '" 

where for the last trace we need in addition fi > 1/2 + 3/2p. Note that the 
embeddings 

E M)U (J) x E M , fl (J) [C(J x fi) n C(J; C 1 ^)]^ 1 , E M)h (J) C(J; C 3 (E)) 

require /i > 1/2 + (n + 2)/2p, which is feasible since p > n + 2 by assumption. 
This restriction is needed for the estimation of the nonlinearities. 

The assertions for the linear problem remain valid for such fj,, replacing E(J) 
by E /J (J), F(J) by F M (J), for initial data subject to (|6.1[) . This relies on the fact 
mentioned above that d/dt admits a bounded iJ°°-calculus with angle tt/2 in the 
spaces L pfl (J;Y). Concerning such time weights, we refer to Priiss, Simonett 
and Zacher [ST] for the Stefan problem, and to Kohne, Priiss and Wilke [TT] 
for the isothermal incompressible two-phase Navier-Stokes problem, and also to 
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Meyries and Schnaubelt [15] for a general theory. Thus as a consequence of these 
considerations we have the following result. 

Corollary 6.2. Let p > n + 2. /i e (1/2 + (n + 2)/2p, 1], a > 0, and suppose 
ip G C 3 (0, oo), /i,d£ C 2 (0,oo) are swc/i i/ia£ 

Kj(s) = -sif)"(s) > 0, /ii(s) > 0, dj(s) > 0, s G (0, oo), i = 1,2. 

Lei f2 C M n fee a bounded domain with boundary dQ, G C 3 ~ and suppose Tq <Z fl 
is a closed hypersurface. Assume the regularity conditions 

(u ,6 ) e [^- 2 / p (n\r )nC(fi)] n+ \ r g w£+ 2 "- a /* 

i/ie compatibility conditions 

div uo = in O \ To , 
Uq = d v $o = in dVL, 

M = Po] = 0, P r „ MVu + [Vwo] T )iT„] = on r 0l 

[^{e )\+aH To =0 on r , 

[d a^o] e % 4 ^- 2 - 6 / p (r ), 

as ?/je£Z as i/ie well-posedness condition 

K@o) 7^ on Tj and 8 > on f2. 

T/ien i/ie transformed problem (|4. 1 [) admits a unique solution (u, n, [7r], 0, /z) G 
E /J ((0, a)) /or some a > 0. TTie solution depends continuously on the data. For 
each S G (0, a), i/ie solution belongs to E(<5, a), i.e. regularizes instantly. 



7. Proof of the Main Result 

In this section we prove Theorem 16.11 for given initial data To G Wp 3 ^ p , 
u G Wp~ 2/p (n \ r ) n and O G Wp~ 2/p {n \ r ) satisfying the compatibility 
conditions and the well-posedness condition stated in Theorem 16.11 We suppose 
that ipi G C 3 (0,oo), fii,di G C 2 (0,oo) satisfy 

Ki(s) = -sip"(s) > 0, fXi(s) > 0, di(s) > 0, s G (0, oo), % = 1,2. 

According to the considerations in Section 4, To can be approximated by a real 
analytic hypersurface for any prescribed r\ > in the sense that d//(A/" 2 S, A/" 2 r ) < 
77, and is parameterized by ho G 3 ^ P (S). It is sufficient to prove the local 
well-posedness of the nonlinear problem (|4.1I) . because the notion of L p -solution 
means that the solution of (|1.1|) is obtained as the push-forward of an L p -solution 
of the transformed problem (14.11) . In order to facilitate this task, we rewrite 
(|4.1[) in quasilinear form, dropping the bars and collecting its principal linear part 
on the left hand side. We set as before ^o(x) = (j,(6o(x)), Ko(x) — k(9o(x)), 
d (x) = d(9 (x)), l (x) = l{9 (x)), and h(t,x) = [V/(e As * 7 (9o(x))]. We then have 
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dfU ~ /Uo(a;)Au + \7tt 


= F u (u,ir,9,h) 


in 




divu 


= F d (u,h) 


in 


n\s, 


M 


= 


on 




— P^f/ynfrW?/ 4- [V?/1 T W^1 
UPO \^) \ v u i I v U J J^SJ] 




on 




-2[ W (s)V!i]f E • ^ s + [tt] - crA s /i 


= G n u ° r {u,e,h) 


on 




w, = 9^6* 


= 


on 


an, 


no(x)d t 6 - d {x)A9 


= F e (u,6,h) 


on 




PI 


= 


on 


s, 


h{t,x)6 + crA E /l 


= G e (0,h) 


on 


s, 


ft/i- [do(a:)M/io(a:) 


= G h (u,e,h) 


on 


s, 


u(0) = wo, 0(0) 


= e a 


in 


ft, 


MO) 


= h Q 


in 


E. 



(7.1) 



The nonlinearities are denned by 

F u (u,e,TT,h) = {fi(6)-/j,(6o))Au + Mi (/i)V?r -(«•(/- Mi (ft))V)tt-(J?(ft) • V)u 
+ - Mi(/i))Vfl • ((J - Mi(/i))Vu + [(J - M 1 (fe))Vu] T )) 

- n(6)(M 2 (h) : V 2 )w - /j,(6)(M 3 (h) ■ V)u + ^{9)M A {h) : Vu, 
(«,/») - Mi(/i) : Vu, 

G^>,ft) = Pe[(mW - /i(0 o ))(Vu + I V "] T )^] 

- ft[/i(0)(Vu + [Vu] T )M (ft)V s /j] 

- P e [m(0)(M 1 (/i)Vu + [Mi(/i)Vu] T )(v£ - M (/i)V E /i)] 

+ [ M ((J - MOVu + [(J - M 1 )Vu] T )(^ E - M V E /i) • ^]M(/j)V e /i, 
GZ or (u, h) = 1(^(6) - Mflo))(V« + [V U ]> S • ^] 

- [/i(0)(Vu + [V U ] T )M (ft)V s ft • ve] 

- [/i(0)(Mi(/i)Vu + [Mi(/i)Vu] t )(i/e - M (/i)V E /i) • ^s] 
+ <7(ffr(h)-A E fc), 

F (u, (9, ft) - (k(6» ) - K(6>))<9 t 6> - (d(0) - d(9 Q ))A9 - d{9)M 2 (h) : V 2 

+ d'(9)\(I - M 1 (h))V9\ 2 - d(9)M 3 (h) ■ WO - n{9){R{h) • V)0 

- n{9)u-{I - M 1 {h))V9 

+ /i(0)((I - Afi(fc))Vu + [(J - Mi(/t))Vu] T ) : (J - Mi(/i))Vu, 
G e (0, ft) = Jifl - W(6)] - a(H r (h) - A s h), 
G h (u, 9, h) = \{d{9)/l{9) - d(9o)/l{9o))d v 9\ + u ■ fa - M (/i)V E /i) 

- [(d//)Mi(/i)Vfl • (^ - Mo(ft)V E ft)] - \{d/l)V9 • M (/i)V E /iJ. 
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Here we employed the abbreviations 

M 2 (h) = Mi(/i) + Mj(h) - M^Mjih), 
M 3 (h) = {I-MxQi)) : VMi(fe), 

M 4 (h) = ((I - M 1 (h))V)M 1 (h) - [((I - A/ 1 (ft))V)Af 1 (/ 1 )] T . 

We prove the local well-posedness of (|7.1[) by means of maximal L p -regularity of 
the linear problem (Theorem I5.2[) and the contraction mapping principle. The 
right hand side of problem (|7.1I) consist of either lower order terms, or terms of 
the same order as those appearing on the left hand side but carry factors which 
can be made small by construction. Indeed, we have smallness of h^, V^/io and 
even of Vf^/io uniformly on E, because Tq is approximated by E in the second 
order bundle. 8 appears nonlinearly in ijj, k, fi, d, but only to order zero; hence e.g. 
the difference fi(6(t)) — /i(#o) will be uniformly small for small times. 

We introduce appropriate function spaces. Let J = [0, a]. The solution spaces 
are defined by 

Ei (a) := {u e H^{J;L p (n)) n n L p {J-,H^{n\T,)) n : u = on dfl, [u] =0}, 
E 2 (a) :=L p {J;H${Sl\'£)), 

Eg (a) := W^-^iJ; L„(E)) D L p (J; W^'p)), 

E 4 (a) := {0 e Hl(J;L p {n))nL p (J;H^{n\Y,)) : d v Q = on 90, [0] = 0}, 
E 5 (a) := Iff ^^(J;!^)) n ^^/^(J;^^)) nL p (J;^-^(E)). 
We abbreviate 

E(a) := {(u, 7T, [tt], 6», /i) e E 1 (a) x E 2 (a) x E 3 (a) x E 4 (a) x E 5 (a)}, 

and equip Ej(a) (j = 1, . . . , 5) with their natural norms, which turn E(a) into a 
Banach space. A left subscript always means that the time trace at t — of the 
function in question is zero whenever it exists. 
The data spaces are defined by 



Fi(a) 


:-L p (J;L p (r!)) n , 


F 2 (a) 


:= H l p {J;H- l {n))nL v {J;Hl{n \ £)), 


F 3 (o) 


:= W p 1 / a -Vap( i7; Lp ( E ))« n L P (J; w£ _1 / p (E)) n 


F 4 (a) 


:=L p (J;L p (n)), 


F 5 (a) 


:= W / p 1 ~ 1 / 2p (J; L P (E)) Pi L p ( J; ^"^(E)), 


F 8 (a) 


:= W p 1 /2-V2P(j. ip(E)) n Lp[J . W^-Vp( S )). 



We abbreviate 

6 

F(a) := {(fu,fd,9u,9e,9e,9h) G J|F 3 (a)}, 

3=1 
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and equip ¥j(a) (j = 1, . . . , 6) with their natural norms, which turn F(a) into a 
Banach space. 

Step 1. In order to economize our notation, we set z — (u,tt, [ivj,0,h) G E(a) 
and reformulate the quasilinear problem (|7.1I) as 

Lz = N(z) (u(Q),6(Q),h(p)) = (<uo,e ,ho), (7.2) 

where L denotes the linear operator on the left hand side of (|7.1|) . and N denotes 
the nonlinear mapping on the right-hand side of (|7.1j) . From Section 5 we know 
that L : E(a) — > F(a) is bounded and linear, and that L : oE(a) — > o¥(a) is an 
isomorphism for each a > 0, with norm independent of < a < ao < oo. 
Concerning the nonlinearity N, we have the following result. 

Proposition 7.1. Suppose p > n + 2, o~ > 0, and let ipi G C 3 (0,oo), fii,di G 
C 2 (0, oo) such that 

Ki(s) = sij}-(s) > 0, fii(s) > 0, di(s) > 0, s G (0, oo), i = 1,2. 

Tften /or eacft a > £fte nonlinearity satisfies N G C 1 (E(a), F(a)) and zis Frechet 
derivative N' satisfies in addition N'(u,ir, [7r|,0, ft) € £(oE(a), oF(a)). Moreover, 
there is r\ > sweft £fta£ /or a given z* G E(ao) ura'ift |fto|c 2 (£) < tftere are 
continuous functions a(r) > and /3(a) > wztft a(0) = /3(0) = 0, sweft iftai 

||iV'(z + z*)|s( E( a ), F(a)) < a(r) + /3(a), I £ B r C E(a). 

Proo/. The nonlinear terms F u (u, 9, tt, ft), F d (w, ft), G* t a "(u, 0, ft) and G£ or (u, M 
are essentially the same as those for the isothermal two-phase Navier-Stokes prob- 
lem (cf. Proposition 4.2 in [H]). The nonlinearities Fg(u,9,h), Gg(0,h), and 
Gh(u,6,h) are similar to those for u = which have been discussed in [21]. In 
fact, all appearing extra terms are of lower order. □ 

Step 2. We reduce the problem to initial values and resolve the compatibilities 
as follows. Thanks to Proposition 4.1 in [11], we find extensions /J G F2(a), 
<?* G F 3 (a) which satisfy 

/J(0) = divu , P s gZ = -P^o(Vu Q + [V^ofKJ- 

Then we define 

g* g := e Ast [G e (0,ft)(O)], := e A ^G h (u ,9 ,h ), 

and set /* = fg = 0. With these extensions, by Theorem 5.2 we may solve 
the linear problem (|5.1|) - (|5.3[) with initial data (ito,#o,fto) an d inhomogeneities 
(/«> /rf) ffui fe^9e' 9h)i which satisfy the required regularity conditions and, by con- 
struction, the compatibility conditions, to obtain a unique solution 

z* = (u*,7r*,[7r*],6>*,ft*) £ E(J) 

with u*(Q) = u 0) 6>*(0) = 6» , and ft*(0) = ft . 

Step 3. We rewrite problem (|7.2[) as 

Lz = N{z + z*) - Lz* =: K(z), z G oE(a). 
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The solution is given by the fixed point problem z = L 1 K(z), since Theorem l5.2l 
implies that L : oE(a) — > o^(a) is an isomorphism with 

\L U( F(o), E(o)) < M, a e (0, do], 

where M is independent of a < ao- We may assume that M > 1. Thanks to 
Proposition 17.11 and due to K(0) = N(z*) — Lz*, we may choose a € (0, ao] and 
r > sufficiently small such that 

r 1 

\K(0) | F(o) < — , \K'{z)\ c{aE(a)tana)) < —, z€ E(a), \z\ E(a) < r 

hence 

\K{z)\ ¥(a) < -L, 

which ensures that L~ 1 K(z) : B oK ( a ) (0, r) — > B oE ( a )(0, r) is a contraction. Thus 
we may employ the contraction mapping principle to obtain a unique solution on 
a time interval [0,a], which completes the proof of Theorem 16. II 

8. The Local Semiflow 

We follow here the approach introduced in Kohne, Priiss and Wilke [TT] for the 
isothermal incompressible two-phase Navier-Stokes problem without phase transi- 
tions and in Priiss, Simonett and Zacher [21] for the Stefan problem with surface 
tension. 

Recall that the closed C 2 -hypersurfaces contained in ft form a C 2 -manifold, 
which we denote by MH 2 (fl). The charts are the parameterizations over a given 
hypersurface S according to Section 4, and the tangent space T-^Ad'H 2 (fi) consists 
of the normal vector fields on E. We define a metric on MW 2 ^) by means of 

<W(£i,£2) :-rf ff (AA 2 S 1 ,AA 2 S 2 ), 

where du denotes the Hausdorff metric on the compact subsets of M". This way 
A41L 2 (n) becomes a Banach manifold of class C 2 . 

Let d,Y,(x) denote the signed distance for S as in Section 4. We may then define 
the level function ip^ by means of 

where 

(j)(s) = s X (s/a) + (sgns)(l - x(s/a)), s e R. 

It is easy to see that S = (^^(O), and Vfs:(x) = vs(x), for each ieS. Moreover, 
k = is an eigenvalue of V 2 (^s(a;), the remaining eigenvalues of V 2 <^s(a:) are the 
principal curvatures Kj of £ at x 6 E. 

If we consider the subset _M% 2 (r2,r) of AAT-L 2 {n) which consists of all closed 
hypersurfaces L G J\4Ti 2 {Vl) such that L C Q satisfies the ball condition with fixed 
radius r > (see the definition below) then the map $ : MU 2 {£l,r) C 2 (n) 
defined by <&(T) — ipr is an isomorphism of the metric space A4H 2 (Q,r) onto 
$(MH 2 (n,r)) cC 2 (ti). 
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Let s-(n-l)/p > 2; for T £ MH 2 (n,r), wedefineT £ W p s (ft,r) fi> r £ Wp S (ft). 
In this case the local charts for T can be chosen of class W£ as well. A subset 
A C r) is (relatively) compact, if and only if C W^(ft) is (relatively) 

compact. 

As an ambient space for the state manifold SM of problem (|1.1[) we consider 
the product space C(ft)™ +1 x .M?{ 2 (ft), due to continuity of velocity, temperature 
and curvature. 

We define the state manifold SM as follows. 

SM :={(«, 0,r) e C(ft)" +1 x MW 2 : (u,9) £ W*- 2/p (Q\T) n+1 , T £ T^ p 4 ~ 3/p , 
div u = in ft \ T, > in ft, u = d v 6 = on 5ft, 
ty(0)] + CTi? r = Prb40)(Vu + [Vu]~> r ] = on T, 1(0) ? on T, 

[d(0)3„0] e v^ 2 - 6 /p(r)}, 

Charts for these manifolds are obtained by the charts induced by .M"H 2 (ft), fol- 
lowed by a Hanzawa transformation as defined in Section 4. 

Applying Theorem 16.11 and re-parameterizing the interface repeatedly, we see 
that yields a local semiflow on SM. 

Theorem 8.1. Letp > n+2, a > 0, and suppose ipi 6 C 3 (0,oo), fii,di £ C 2 (0, oo) 
such that 

m(s) — -sip"(s) > 0, fii(s) > 0, <2j(s) > 0, s € (0, oo), i = 1, 2. 

TTien problem generates a local semiflow on the state manifold SM. Each 

solution (u,9,T) exists on a maximal time interval [0,£*), wheret* = £*(«o, To). 

Note that the pressure does not occur explicitly as a variable in the local semiflow, 
as the latter is only formulated in terms of the temperature 9, the velocity field 
u, and the free boundary T. The pressure tt is determined for each t from (u, 9, T) 
by means of the weak transmission problem 

(V7r|V^) i2( n) = 2(div(/z(0) J D) - u • Vu|V^) Za(n ), <P G i^(ft), 

W = 2[ M (0)D] on T. 

Concerning such transmission problems we refer to Section 8 in 

Let (it, 0, r) be a solution in the state manifold SM with maximal interval of 
existence [0,i*). By the uniform ball condition we mean the existence of a radius 
ro > such that for each t £ [0,£«), at each point x £ T(t) there exists centers 
Xi £ fti(£) such that B ro (xi) C fti and T(i) n -8^,(2^) = {x}, i — 1,2. Note that 
this condition bounds the curvature of T(t), prevents parts of it to shrink to points, 
to touch the outer boundary 9ft, and to undergo topological changes. 

With this property, combining the local semiflow for ([Lip with the Ljapunov 
functional and compactness we obtain the following result. 
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Theorem 8.2. Letp > n+2, a > 0, and suppose ipi <E C 3 (0,oo), fii,di G C 2 (0, oo) 
such that 

Ki{s) = -sip"(s) > 0, m(s) > 0, di(s) > 0, s e (0, oo), i = 1,2. 

Suppose that (it, 9, T) is a solution of in the state manifold SM. on its maxi- 

mal time interval [0, t*). Assume there is constant M > such that the following 
conditions hold on [0,£*): 

(i) \u(t)\ W 2- Vpn , \e(t)\ W 2-2 /P , \r(t)\ W 4- 3/P , \\d{e{t))d v u(t)\\ W 2- e/P < m«x> ; 

(ii) \l(e(t))\,6(t) > l/M; 

(iii) T(t) satisfies the uniform ball condition. 

Then t„ = oo, i.e. the solution exists globally, and uj + (u,9,T) C £ is non-empty. 

Proof. Assume that (i), (ii), and (iii) are valid. Then r([0,i*)) C Wp~ 3/p (n, r) 

is bounded, hence relatively compact in Wp e (0, r). Thus we may cover 
this set by finitely many balls with centers real analytic in such a way that 
dist^-s/p-* (r (t),T,j) < 5 for some j = j(t), t £ [0,t*). Let J k = {t e [Q,t*) : 

j(t) = fc}; using for each k a Hanzawa-transformation Hfc, we see that the pull 
backs {(u(t, -),6{t, ■)) o S fc : te J k } are bounded in W / p 2_2/p (rj \ S fc )™ +1 , hence 
relatively compact in Wp 2 ^ p £ (Q \ Sfc)™ +1 - Employing now Corollary 16.21 we 
obtain solutions (it 1 , 6 1 , T 1 ) with initial configurations (u(t), 9(t), T(t)) in the state 
manifold on a common time interval, say (0, t], and by uniqueness we have 

( U x (r), x (t), r x (a)) = (u(t + r), 0(t + r), r(* + r)). 

Continuous dependence implies that the orbit of the solution (u(-),6(-),T(-)) is 
relative compact in SM, in particular t* — oo and (m, 0, r)(M+) C SA4 is relatively 
compact. The negative total entropy is a strict Ljapunov functional, hence the 
limit set u>+(u, 9, T) C SM. of a solution is contained in the set 8 of equilibria. By 
compactness w+(it, 0, T) C SM. is non-empty, hence the solution comes close to £ , 
and stays there. □ 

Remark. We can prove that each equilibrium is normally hyperbolic. Therefore, 
taking into account the results in |21j . one may expect that each orbit satisfying 
the assumptions of Theorem 18.21 actually converges to an equilibrium. This topic 
will be considered elsewhere. 
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